Multi-optimization for thermal deformation of gravitational wave telescope based on CFRP characteristics

Gravitational wave telescope place extremely high demands on structural thermal deformation, making material selection a critical issue. Carbon fiber reinforced polymer (CFRP) is an ideal choice for the support structure of telescope due to its low coefficient of thermal expansion (CTE) and designable properties. However, current research on the optimization of the CTE of CFRP is scarce, and conventional methods struggle to find layups that meet the requirements. In this paper, an unconventional layup optimization method is proposed to solve this problem. Initially defining the characteristics of the telescope structure and using different layup material for the main and side support rods to minimize thermal deformation. Subsequently, the NSGA-II algorithm is used to optimize the layups which are divided into conventional and unconventional layups. Specimens are then produced from these results and tested to assess the impact of processing errors on practical applications. The results demonstrate that the optimized CFRP meet the CTE requirements and, when applied to the structure, significantly reduces the thermal deformation in the eccentric direction compared to conventional designs. Additionally, a numerical analysis evaluates the effect of ply orientation errors on the performance of unconventional layups, discussing the method's limitations within these contexts.

There are geometric relationships: The basic equation for thermal expansion is �L = L × α × �T , which is inserted into the Eqs.(1) and (2)  to obtain the thermal deformation of the node U as a function of the geometric parameters of the structure: In the general design of truss structures, the same material is typically chosen for all support rods, denoted as α s = α P .However, it has been found that the structural thermal deformation, dy, can be reduced by assigning different values such that α S = α P , as demonstrated by Eq. (3).The principle of the optimized design is shown in Fig. 3.
Based on the current design of the telescope prototype, the basic structural design parameters are as follows: L A = 326 mm, L B = 215 mm, L C = 838 mm, and the global space temperature change of the telescope during the orbital phase is ΔT = 1 K.
The objective of the structural design is to minimize the thermal deformation of the structure, and the observation of Eq. (3) reveals that changing the thermal deformation dy in the direction of deviation is achievable with different parameter pairings to achieve dy = 0. Therefore, the parameters with dy = 0 are inserted into Eq.( 3) to obtain the short form: (1) �L S = �L C cos β × cos γ (2)  The laminates' in-plane x-direction CTE, representing the axial thermal expansion of the rods, is denoted by α x , while the radial thermal expansion is represented by α y .Notably, the joints in the structure are all made of Invar, with CTE of α Invar = 1 × 10 −6 /K .Thus, a third design indicator emerges: it is crucial not only to ensure that α x in each support rod meets the specified requirements, but also to adjust α y to be as close as possible to 1 × 10 −6 /K .Figure 4 illustrates the orientation correspondence between the support rods and the material.The simulation in this figure uses millimeter units to analyze the thermal expansion of CFRP laminates and support rods per unit length (L = 1 m) under conditions of unit temperature ( T = 1 K ) change.
In summary, given the need to perform layup optimization design separately for the main and side rods, there are four optimization objectives: α S x , α S y , α P x , α P y .Thus, there are three requirements for the CTE of materials: 1) α S x : α P x = 1:1.217;2) α P x < 0.6 × 10 −6 /K ; 3) α P y , α S y as close to 1 × 10 −6 /K.www.nature.com/scientificreports/

Mathematical modeling of thermal deformation of CFRP
From the above analysis, it is clear that the core of the thermal deformation of the structure lies in the CTE of the three truss rods.Since this value can be altered by designing the angle of the CFRP material layup, it is necessary to establish a functional relationship between the CTE of laminates and the orientation of the lamina.CFRP is obtained from a single layer of plate after lamination.The laminates can be considered as a beam structure with continuous variation in the thickness direction, as shown in Fig. 5.
Based on the diagram, it can be known: Strain formula can be obtained: From the structural strain, the structural stress of each layer can be deduced in reverse, Q k is the off-axis stiffness matrix for each layer: By integrating the stress of each layer in the thickness direction, the following can be obtained: where N x and M x are the resultant force and moment per unit length on the laminated composite material in the x direction of the plane.The same applies to the y direction and the shear direction x-y.
Obtaining the constitutive relationship of laminated composite materials: A ij is the stiffness coefficient that relates the resultant force to the mid-surface strain; D ij and B ij indicate the coupling of the orthogonal deformation to other directions.
The objective of this work is to optimize for the thermal deformation of the material.The reason is that the thermal expansion does not involve bending stresses and the order of stacking of constraint layers is symmetrical, which prevents the coupled deformation of the material tension-bending-shear, so that B ij = D ij = 0 .Therefore, Eq. ( 10) is simplified and inverted:  A ij indicates that the element is from the inverse matrix of the original tensile stiffness matrix [A] ; N T rep- resents the thermal stress per unit width experienced.Since the model of the thermal expansion of the material does not take external forces into account, the stress of each layer is The relationship between the N T and CTE of each layer is determined as follows: Meanwhile, the relationship between the elements in the tensile stiffness matrix [A] and the stiffness matrix of each layer is: The material strain vector [ε] for a certain temperature difference T is obtained by substituting [N T ] and A ij in Eq. (11), and the relationship between the coefficient of thermal expansion of the laminate [α] l and the orientation of the lamina is obtained by integration: [t] is the coordinate transformation matrix.Due to the symmetry of the layup, thermal expansion occurs only in the orthogonal directions, and the coupling thermal expansion coefficient α xy = 0 .Therefore, Eq. (14) shows that α x and α y within the laminate are adjusted by the orientation of the lamina.

Optimizing CFRP layup by NSGA-II
This paper performs layup optimization in two stages: first, considering the lower CTE of the side support rods than the main rod, the layup for the side support rods is optimized to achieve the minimum α S x , α S y ; then, the layup of the main support rod is optimized in a ratio of 1:1.217 to obtain α P x , α P y , as shown in Fig. 6.To meet the requirements for structural strength and geometric dimensions, 20 layers of 0.15 mm thick M40JB/C1414 CFRP are used.Table 1 shows the performance of the alternative materials for the telescope structure design.The Zerodur was supplied by SCHOTT with the specification of ZERODUR® Expansion Class 1. Invar alloys was supplied by Beijing Aerospace Invar Materials Technology Development Corporation, which is a kind of modified alloy based on 4J32.CFRP was supplied by Chang-Chun Long Aerospace Composite materials Co., Ltd. with the fiber specification M40J, and the matrix material was selected as C1414 epoxy resin.
To prevent twisting and bending coupling effects, the laminate should be laid in a complementary and symmetrical manner.The first 10 layers and the last 10 layers must be symmetric.Within these 10 layers, it is necessary to ensure the presence of layups with identical positive and negative angles.This arrangement requires optimization for five ply angles, referred to as (±θ 1 , ±θ 2 , ±θ 3 , ±θ 4 , ±θ 5 ) s , where θ i ∈ −90 • , 90 . Layer optimization process.
Table 1.Candidate materials for truss support structure of telescope prototype.www.nature.com/scientificreports/optimization of layup design involves directly listing several common layups and selecting the one that is closest to the design goal as the design scheme.This method is simple and convenient but lacks precision.Therefore, this paper presents optimization of conventional and unconventional ply angles.Conventional ply angles are fixed at [± 45°, ± 30°, ± 60°, 0°, 90°] and can be controlled by variable constraints in the optimization program.Unconventional layups are not constrained during optimization.
The study utilizes the NSGA-II algorithm, which requires conflicts between different optimization objectives to find a set of design solutions forming a 'Pareto front' .The results on the Pareto front represent extreme states, where improving one objective function often implies a regression in another, ensuring that all optimization results are relatively optimal solutions.For laminate design, α x and α y are conflicting objective functions, making the NSGA-II algorithm particularly suitable for optimizing the CTE issue of CFRP.
Moreover, the advancement of NSGA-II lies in its introduction of the concepts of ' crowding distance' and 'non-dominated sorting' , ensuring the effectiveness and diversity of the optimization results.
This paper presents a program for optimizing CFRP layup using Python.The program consists of 6 optimization steps, as shown in Fig. 7.

Initial setup and problem definition:
The initial population for the optimization scheme is composed of 100 randomly generated layup design schemes, each referred to as an individual.These individuals contain five optimization variables, referred to as genes, which are ply orientations (θ i , i = 1,2, 3,4, 5) .The subsequent selection, crossover, and mutation processes are carried out using these genes.The termination criterion for the process is set to 120 iterations to achieve the Pareto front.2. Nondominated sorting categorizes the population into different non-dominated layers based on the dominance relationship between individuals.For example, if individual P has smaller values of both α x and α y compared to individual Q, then Q is considered to be dominated by P. The number of individuals, including P, that can dominate Q within the population is tallied, denoted as f (Q) = n .In short, genetic prioritization between individuals is determined by comparing their performance across individuals.In this study, the performance manifestation is the CTE α x and α y of the different layup schemes.
As shown in Fig. 8a: For individual A, there are no other individuals whose α x and α y are both lower than A's, meaning there is no individuals that can dominate A, resulting in f (A) = 0 .Conversely, individual F is dominated by 2 individuals, B and C, whose α x and α y are both lower, denoted as f (F) = 2 .However, observing the layup scheme F, compared to layup schemes A and D, only one of α x or α y performs better, neither of which has a dominance relation.Similar assessments apply to the rest of the population.
Individuals with the same dominance counts are grouped into equivalent sets and ordered by their dominance count.These sets are labeled sequentially as the 1st Pareto front, 2nd Pareto front, and so forth, until all non-dominated Pareto fronts within the population are delineated.This procedure identifies the most preferable layup designs, positioning those within the earliest Pareto fronts as high-priority design options, earmarked for preferential consideration in subsequent genetic operations.As illustrated in the figure, D(B) > D(A), indicating that although the layup schemes A and B have the same non-dominance ranking and are both on the 1st Pareto front, the priority of layup B is higher than that of A. This ensures a relatively uniform distribution of results within a certain range.4. Selection, crossover, and mutation operations of genes: Layup design schemes with high comprehensive priority are selected through the non-dominated sorting and crowding distance comparison.Subsequently, crossover and mutation operations are performed on their five ply orientations (θ). 5. Merging the population and optimizing iteratively: After performing selection, crossover, and mutation, new layup designs are obtained.These new and existing layup schemes are then merged, and non-dominated sorting and crowding distance calculations are performed again.This process is iterated until the optimal Pareto front is reached, as shown in Fig. 9. 6.The obtained output results yield several layup designs, all of which are relatively optimal.Only those layups that meet the requirements for α x and α y should be incorporated into the structural design.

Experiment and result
Considering that each ply orientation may deviate from the design value in practical applications, potentially affecting material performance, it is necessary to conduct validation experiments.These experiments should analyze errors in layup angles in conjunction with actual manufacturing processes.This approach will facilitate a more in-depth discussion of the applicability of the CFRP optimization method in ultra-stable structure applications.

Optimization design results
Select the layup from the optimization results conducted using the NSGA-II algorithm, as shown in Table 2.The first four layups are conventional angle layups, which struggle to meet all three design objectives simultaneously due to the angle constraints.Layups 5-8 feature unconventional layups that are not constrained by angles during the optimization process, thus allowing for a broad range of feasible design results that meet the requirements.The last two layups, designed to demonstrate the superiority of the algorithm, achieve extremely low axial CTE, ( 15)  www.nature.com/scientificreports/illustrating the method's effectiveness in optimizing CFRP.However, due to their excessive lateral thermal expansion, these layups are not practical for structural applications and are therefore not discussed in detail.Layup 8 has been identified as meeting the specified requirements for both axial and radial thermal expansions more effectively.This design ensures that the axial CTE is less than 0.6 × 10 −6 /K and maintains a lower radial CTE of α S y = 2.736 × 10 −6 /K , which facilitates a better thermal stress match with Invar joints.When selecting layup 8 for the laminates used in the side support rods, to meet the requirement of α S x : α P x = 1 : 1.217 , the CTE of the main support rod should be α P x = 0.4322 × 10 −6 /K .The results of the second stage optimization specify layup 7.

Experimental results and error analysis
Given that angle errors in lamination can affect material properties, the actual structural performance might deviate from the intended design outcomes.To ensure the optimized layups perform adequately in applications, sample tests were conducted on the first eight layups to evaluate their CTE.
Considering the lack of precision in roll-wrapped CFRP prepreg for unconventional layups due to non-integer layup angles, choosing the winding process for manufacturing is more appropriate.This winding process, in which the angle of layup is controlled by controlling the rolling and traveling speed of the equipment, not only reduces errors, but also ensures uniform material properties, which is particularly effective for processing CFRP pipe fittings.The CFRP tubes were then cured, demolded, and post-treated, and their ends were cut to standard dimensions based on the method in the references 21 to obtain the CTE test specimens.The test sample production flow is shown in Fig. 10.
As can be seen from Eq. ( 14) in the sectoin "Mathematical modeling of thermal deformation of CFRP", it is important to obtain basic lamina mechanical and thermal expansion data before proceeding with the layup design.Therefore, mechanical property parameter tests and CTE tests were conducted for the selected M40J/C1414 material.The Fig. 11a is the sample of tensile test of lamina material based on the standard GB/T 3354-2014.The Fig. 11b is the CTE test of lamina material and 8 types of layups based on the testing standard GB/T16535-2008.Considering the anisotropic characteristics of CFRP, all of the above include test results in both longitudinal and transverse directions.
The CTE of specimens from three different batches was measured through experiments using the non-contact LINSEIS L75 Laser Dilatometer.Each sample underwent three measurement cycles at temperature conditions ranging from 20 to 80 °C.The average CTE values, calculated from these three cycles, are presented in Tables 3  and 4.
Here is the error and uncertainty analysis of 8 types of layups, as shown in Fig. 12.
The test results of CFRP specimens with common ply orientations show more stability, indicating that the current level of lamination technology for conventional angles such as ± 45°, ± 60°, ± 30°, 0°, and 90° is superior to that of unconventional angle layup designs.In terms of thermal expansion performance, the CTE is smaller for specially designed unconventional layups.However, these unconventional layups schemes are more challenging to achieve when considering manufacturing errors.

Application in gravitational wave telescope design
Before conducting thermal deformation analysis of the structure, data of CFRP from 1st batch of layups 7 and 3rd batch of layup 8 was integrated into the structure.As demonstrated in Fig. 13, a comprehensive finite element simulation was conducted using ABAQUS.The general performance results are presented in Table 5.  Solid element (C3D8R) meshes were used for the entire structure, excluding the three CFRP support rods.The support rods are modeled as an anisotropic material shell element (S4R), with the basic mechanical properties of the material consistent with theoretical values, and the CTE chosen based on actual test results.The material definition for the CFRP support rods is presented in Table 6.
In order to accurately obtain the deformation of the truss, the structure support metal frame was constrained to prevent any participation in the deformation.Thermal deformation was evaluated based on the displacement data at Node U, located at the secondary mirror's position, as detailed in Table 7.
Due to the high specific stiffness of the CFRP material and the overall light weighting of the structure, the three items of structural weight, modal and deformation under gravity are satisfied with the design requirements.
It should be noted here that the maximum stress indicator presented in row 4 is based on the strength requirements of the M40J/C1414 lamina material.Given that the reinforcement of fibres in the transverse direction is relatively modest, the strength of CFRP material is the lowest in that direction, with a value of 45 MPa.Given the strength margin, it is deemed appropriate to set the maximum strength design target at 35 MPa.A finite element simulation analysis was conducted to assess the working condition of the telescope structure prototype.The analysis revealed that when thermal stress is combined with the influence of gravity, the maximum stress of the structure as a whole is 4.96 MPa.Consequently, the current structural design is deemed to meet the maximum strength requirements.
It is clear that the performance of CFRP in actual applications will deviate from the design values, and such deviations introduced into the structure can have unexpected effects on the thermal deformation of the telescope.Therefore, it is necessary to incorporate the actual measurement results of CFRP into the telescope structure, compare them with the theoretical design results, and determine the acceptable range of errors.
Table 7 shows the optimization results for thermal deformation of the telescope structure.The material of the three rods was replaced to compare different design schemes, while keeping all other components unchanged.
The optimized structure achieved the best performance under ideal conditions, exhibiting the lowest thermal deformation compared to all other scenarios.It reduced dz by 14% and dy by 98% compared to the all-Invar  structure.Furthermore, an examination of rows 2 and 3 reveals that the optimization of thermal deformation in the structure through the use of CFRP, with a CTE ratio of 1.217:1 for the main support rods to the side support rods, is evident.This is in contrast to the unoptimized quasi-isotropic layup, where dz was reduced by 62.05% and dy by 99%.Observing rows 3 and 5, introducing materials with error into the structural simulation will indeed affect the thermal deformation performance and make the structure fail to achieve the desired thermal deformation results.However, in comparison, the optimization effect is still significant and the structure still meets the design requirements.

Discussion of limitations of optimization methods
Following the analysis presented, it is apparent that deviations in the CFRP layup have a substantial impact on the material's actual performance.Therefore, an in-depth analysis should be conducted in conjunction with the error numerical analysis method to identify the limitations of the optimization method proposed in this paper.
Given the significant impact of the axial CTE of the main support rod on the structure's thermal deformation, the discussion will henceforth focus solely on α P x , hereinafter referred to as α P , to control variables.Using the Design of Experiments (DOE) module within Isight, these stochastic errors were integrated into the previously determined layup of the main support rod, providing 360 random error calculations and statistically determining the possible outcome range for α P .The examination of the implications of a 5° lamination error is illustrated in Fig. 14, with similar analyses applicable to other error scenarios.
When the error is within ± 5°, the resulting α P (10 −6 /K) interval is [-0.31, 1.12].The theoretical prediction for α P is 0.434 as layup 7 in Table 3, and the experimental test results for three different batches of samples are 0.83, 0.98, and 0.51, all within this range.This indicates that the experimental results align with the error prediction, albeit with significant dispersion.To further investigate the applicability of the optimization method presented in this study, the results of the DOE under the conditions of each ply orientation error of ± 2.5°, ± 1°, ± 0.5°, and ± 0.1° are presented, as shown in Fig. 15.The higher 3σ boundary values from the DOE numerical experiments are incorporated into the truss model to analyze and obtain dy values for different error conditions, as shown in Table 8.
The results show that:    1.When the layup error level is controlled within 0.1°, the optimization design method proposed in this paper significantly reduces structural thermal deformation.Although a substantial difference remains compared to the ideal state, the method sufficiently meets the design requirements of dy < 200 nm/K and dz < 550 nm/K, indicating its effectiveness.2. The optimization method is no longer applicable when the error level exceeds 0.5° or above.3. When the error exceeds 1°, there is a likelihood that the structural performance will not meet the design requirements.
It is clear that utilizing layup optimization to adjust α P = α S for minimizing structural thermal deformation demands stringent precision in the CTE of the rods.It is worth noting that even a slight deviation of 0.1° in each layer can result in actual outcomes that deviate significantly from the ideal design solution, by an entire order of magnitude.Currently, based on actual production experience and the results of tests carried out by the LISA team on CFRP in the literature 5 , the typical processing error for unconventional layups is approximately 0.5°.Therefore, to utilize the structural optimization method proposed in this paper to reduce the thermal deformation of the truss structure, the level of lamination technology should be considered, and a higher precision lamination process should be employed in the production of the CFRP truss.

Conclusion
This paper presents three design objectives for the CTE of CFRP, derived from the design requirements of the gravitational wave telescope truss structure.The subsequent sections describe the optimized design method, experimental verification, application in structure, and numerical error analysis.The following conclusions are drawn: 1.This study introduces a novel design method to address the thermal deformation challenges in the structure of gravitational wave detection telescopes, leveraging the properties of CFRP materials and truss structures.
Initially, a mathematical model for the thermal deformation of a three-rod truss structure is established.The method differentiates the CTE between the main and side support rods, using the structural characteristics to mitigate thermal deformation in the eccentric direction of the telescope.Subsequently, considering the structural requirement for a specific ratio of CTE values, the NSGA-II algorithm is employed alongside the design characteristics of CFRP laminates to optimize the axial and radial CTEs, yielding a series of CFRP laminations that meet the design requirements.Afterwards, the design results are divided into conventional and unconventional laminations for sample testing and error analysis, and the most suitable laminations are selected for structural design.The findings confirm that the optimization method proposed in this paper is feasible under certain process conditions and the optimization effect is obvious.2. It is recommended that unconventional angle laminates be employed to address the issue of the material's requirement for a fixed value of CTE.Comparing the test results of conventional layups and unconventional layups, it can be seen that although the processing error of unconventional layups is larger than conventional layups, the thermal expansion performance of the unconventional layups can meet multiple design goals simultaneously.Therefore, the unconventional layups have great development potential for structures that do not require high mechanical efficiency.3.In terms of structural applications, the optimized CFRP design can be used to meet the different CTE requirements of the main and side support rods in the structure, which can significantly reduce the thermal deformation of the structure.This optimization not only enhances the dimensional stability and operational reliability of gravitational wave detection telescopes but also serves as a valuable case study for enhancing www.nature.com/scientificreports/similar structures in aerospace engineering.This is particularly relevant for telescope structures that demand stringent thermal deformation control.4. Considering the large processing error of unconventional laminates in actual situations, this causes the actual CTE of laminates to deviate from the design value, thereby affecting the structural performance.Within the limitations of the study, this paper uses the DOE module in Isight to perform a numerical simulation of the laminate error on the axial thermal expansion coefficient of the selected laminate, with the objective of determining the applicable scope of the method proposed in this paper.5.This paper has some shortcomings and limitations in its research content.It primarily focuses on principle prototype tests of the gravitational wave detection telescope, with experimental objectives centered on detecting thermal deformation of the structure.Consequently, discussions from material optimization to structural design concentrate primarily on thermal deformation, paying less attention to strength, stiffness, light weighting, and other mechanical properties.Additionally, discussions about errors in CFRP layup angles have relied extensively on numerical simulations, as outlined in section "Discussion of limitations of optimization methods", without experimental validation.Future research could benefit from integrating tools like metallographic microscopes and employing techniques such as machine learning and visual imaging to detect errors in CFRP layup angles, thereby enhancing the performance stability of carbon fiber products.

Figure 2 .
Figure 2. Projections of the truss model on the QMU (a) and YOZ (b) planes.

Figure 5 .
Figure 5. Schematic diagram of the deformation of a laminated composite material structure as a beam structure.

Figure 9 .
Figure 9. Iterative optimization to the optimal Pareto front.

Figure 11 .
Figure 11.Layup samples for testing mechanical properties of lamina (a) and CTE tests for laminates (b).

Figure 12 .
Figure 12.Error and Uncertainty Analysis of CTE of Samples in α x and α y .

Figure 13 .
Figure 13.Structural thermal deformation after bringing experimental results into finite element simulation.

26 B 66 B 11 B 12 B 16 B 12 B 22 B 26 B 16 B 26 B 66 D 11 D 12 D 16 D 12 D 22 D 26 D 16 D 26 D 66
12 A 16 A 12 A 22 A 26 A 16 A 26 A 66 B 11 B 12 B 16 B 12 B 22 B 26 B 16 B

Table 3 .
Experimental results of CFRP samples of α x .

Table 4 .
Experimental results of CFRP samples of α y .

Table 5 .
Summary of structural mechanical properties.

Table 6 .
Material definition for the support rods. E

Table 7 .
Results of structure thermal deformation for different cases.